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ABSTRACT 


A quantitative, computer-aided study was made on the spectral separa- 
bility of timber types and condition classes in the Southeastern United 
States, using LANDSAT-1 multispectral scanner data. Conclusions were 
obtained on accuracies at different levels of mapping detail and the choice 
of parameters affecting mapping accuracies, such as spectral bands, number 
of bands, and seasons of data. 

It was concluded that LANDSAT-1 could be used effectively to discrimi- 
nate the gross forest features of softwood, hardwood, and regeneration. 

The only significant detectable age difference would be between an estab- 
lished forest versus a young (or denuded) forest, i.e., regeneration. The 
red or near infrared bands would be better for discrimination; phenological 
early and late spring data would be better than winter (summer and autumn 
data were not available for analysis). Anc a temporal analysis would be 
superior to single-season analysis. Lastly, two spectral bands would be 
most cost effective for computer analysis. 

The study site, Sam Houston National Forest of East Texas, is a typical 
forest in the Flatwoods Zone, Southern Region, U. S. Forest Service. The 
widely accepted computer scheme of training-field, maximum likelihood- 
classifier was employed, while cross -classification accuracies and diver- 
gence measures were computed to evaluate timber type separability. 


1.0 INTRODUCTION 


This study was undertaken to determine the feasibility of mapping 
timber types and condition classes in the Southeastern United States using 
Land Satellite (LANDSAT-1, formerly Earth Resources Technology Satellite, 
ERTS-1) multispectral scanner (MSS) data via computer-aided analysis. 

Timber types refer to stand composition by dominant species, and condition 
classes refer to age, size, health, and adequacy of stocking. To this 
end, three objectives were pursued: 

• To quantitatively determine the spectral separability between timber 
types and condition classes, and additionally between the general 
classes of merchantable timber which consist of sawtimber, pole- 
timber, and regeneration softwood stands. 


*The material of this paper was developed under NASA Contract NAS 
9-12200 and prepared for the Earth Observations Division, NASA, Johnson 
Space Center, Houston, Texas. 
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• To determine the optimal spectral bands and seasons of MSS data 
for maximum separability between timber types and condition classes 

• To determine the effects on classification accuracies because of 
variations in the number of spectral bands 

Such a study is essential to computer-aided remote sensing applica- 
tions on timber resource inventories concerned with mapping and volume 
estimation. The mapping accuracy, mapping levels, sensor spectral band 
coverages, seasons of data acquisition, and number of bands used in analysis 
are important design criteria in these inventories. An optimally designed 
timber resource inventory using remote sensing data, from satellites and 
aircraft, in conjunction with ground surveys would provide a most efficient, 
timely, accurate, and economical solution to many forest management prob- 
lems. For example, studies by^ Aldrich (1971) and Colwell (1973) have proved 
the efficient use of a multistage sampling scheme to estimate timber volume 
on a forest-wide basis. 

The present study determined whether or not satellite MSS data such as 
that of LANDSAT-1 could be used at or beyond the mapping level of softwood, 
hardwood, mixed softwood -hardwood, and regeneration, the four categories 
constituting forest land. (See hierarchy levels in section 2.2.) Also, 
this study established the feasibility of discriminating merchantability 
age classes by LANDSAT-1, i.e., discrimination between sawtimber, pole- 
timber, and regeneration stands of softwood. These conclusions were derived 
from analyzing a typical forest in the Southeastern United States, as 
represented by the Sam Houston National Forest of East Texas. 

A number of past studies have been carried out on similar subjects, but 
not to the same extent of details in mapping and not to the same amount of 
quantification of separability. In the report by Heller et al. (1974), 
investigators reported on the use of October LANDSAT-1 MSS data and computer 
classification techniques in mapping level II land use classes in Georgia. 
Level II included pine and hardwood classes and were identified with 
accuracies ranging from 42 to 81 percent. Erb (1974) reported the analysis 
of August 1972 LANDSAT-1 data on the Sam Houston National Forest of Texas, 
breaking forest land down into hardwood versus pine with accuracies as high 
as 91 percent. Aircraft MSS data have also been analyzed with the purpose 
of identifying forest land use classes. Weber (1972) reported processing 
a November data set over Atlanta, Georgia, and found that spectral bands 
very similar to band 7 (0. 8-1.1 micron) and band 5 (0.6-0. 7 micron) of 
LANDSAT-1 were the third and fourth best channels after the infrared band 
(1.5-1. 8 micron) and thermal band (9.3-11.5 micron) for separability of 
forest classes. Accuracies for classification of pine and hardwoods ranged 
from 20 to 80 percent. 

The present study is part of the Forestry Applications Exploratory 
Studies Project, which is conducted by the Earth Observations Division at 
the Lyndon B. Johnson Space Center of the National Aeronautics and Space 
Administration and by the Southern Region of the Forest Service, U. S. 
Department of Agriculture. Project details can be found in Anon. (1974). 
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2.0 STUDY SITE AND FEATURES 


2.1 Site Description 


The study site is the Conroe Unit, Raven District of Sam Houston 
National Forest, located 90 kilometers north of Houston, Texas (figure 1). 
This forest is in the "East Texas Piney Woods" or "Flatwoods," which is 
the heavily forested portion of East Texas. The Conroe Unit is within the 
southern Gulf Coastal Plain, and the overall area slopes to the southeast 
at 1.5 meters per kilometer. Slopes are generally between 3 and 7 percent; 
elevation differences between stream bottoms and ridge tops are usually no 
more than 25 meters. 

The 28,500 hectares of the Conroe Unit consist of approximately 10 
percent shortleaf pine ( Pinua echinata Mill.) on the ridges and upper 
slopes, 75 percent loblolly pine (Pinus tueda L.) and hardwoods on the 
lower slopes, 10 percent hardwoods in the drainage ways, and 4 percent 
regeneration areas. The remaining 1 percent is made up of openings such 
as pipelines or oil well sites. The most common hardwood types are mixed 
oaks — laurel oak ( Quercus lauri folia Michx.) and willow oak (Queraua 
phellos L.), and gums and oaks - sweetgum ( Liquidambar styraciflua L.), 
nuttal oak ( Quercus nutalli Palmer), and willow oak. 


2.2 Analysis Levels 


A five- level hierarchy of land features was used for this study 
(table I). The terminology and structure of the hierarchy resulted from 
modifications of those from the Society of American Foresters (Ford- 
Robertson, 1971; Anon., 1974), The Forest Survey (Sternitzke, 1967), and 
the Geological Survey (Anderson et al., 1972). The definitions of levels 
I and II features are adopted as follows: 


Forestland 


Softwood - 


Hardwood - 


land of 0.4 hectares or more in size supporting 
a stand of trees whose crowns cover more than 
10 percent of the area 

gymnosperms, generally having evergreen and 
needle foliage; a softwood stard consists of 
more than 50 percent pint in the overstorv. 

angiosperms, generally having broadleaved and 
deciduous foliage; a hardwood stand consists of 
less than 25 percent pine in the overstory 


Mixed softwood-hardwood - a stand of mixed softwood-hardwood consisting of 

25 to SO percent pine and the rest hardwood 

Regeneration area - cutover, burned, or otherwise denuded forestland 

in process of being reclaimed by a young forest. 

The class "regeneration" in level II was considered appropriate in this 
computer-aided remote sensing application, since regeneration has spectral 
characteristics known to be distinct from established forest stands. Also, 
it was considered impossible to discriminate among shortleaf regeneration, 
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loblolly regeneration, and hardwood regeneration. The breakdown under 
regeneration into seedling-sapling or nonstocked is at level V. 

Levels III, IV, and V are traditional breakdown and represent levels 
of detail usually mapped by ground survey or photointerpretation. The 
special level "general age class" was created to study the condition 
classes of softwoods collectively, i.e., sawtimber, poletimber, and regen- 
eration areas. The identification of these general classes, as opposed to 
the more detailed level IV classes, has mercantile and management import. 
The limited amount of hardwood sites in the study area did not allow this 
study to analyze hardwood sawtimber, poletimber, and regeneration separa- 
bility. 


3.0 TECHNICAL APPROACH 


The overall approach taken in this study is the "training-field 
classification/analysis" process used in computer systems such as LARSYS 
(Phillips, 1973), RECOG (Ells et al., 1972), and ERIPS (Anon., 1973). This 
approach essentially consists of (1) the acquisition of spectral signatures 
of land features by locating "training fields" of these features on MSS 
data and computation of their statistics; (2) searching for optimal sec 
of channels for classification through calculation of some mathematical 
distances, e.g., divergence (Marill and Green, 1963) or probability of 
correct classification (PCC) (Anderson, 1958); and (3) computer classifica- 
tion, e.g., using a maximum likelihood classifier which assigns to an un- 
known data pixel (picture element in MSS data) a most likely class associa- 
tion from all the possible training classes. A training class means land 
feature on which training field statistics are available. Through this 
process, timber type separability information will be obtained via the two 
mathematical measures of divergence and PCC. (See section 3.2 for details.) 

Figure 2 is a schematic flow diagram of the analysis procedures 
followed in this study. Three LANDSAT-1 data sets were checked for data 
quality before they were registered and composed to form one 11-channel 
data set. The registration was performed image-to- image and to the 
corresponding longitude and latitude locations. Since each individual - 
date data set has four MSS channels, the composition of data sets resulted 
in a 12 channel data set with the ordering of the channels and their 
spectral coverages tabulated in table II; however, the last channel of the 
May data set was too noisy, and thus omitted from analysis, lienee an 11- 
channel data set. 

Random training fields of the level V forest features (table I) 
were selected and their coordinates recorded on the MSS data. The 
same locations applied to all three dates because the three data sets had 
been registered to one another. By straightforward aggregation, training 
fields for all hierarchy levels were compiled. At each level, divergence 
calculations were made, compiled, tabulated, and analyzed. Also class 
pairs were classified, a pair at a time, producing pairwise correct classi- 
fication accuracies. The average of these were then calculated, tabulated, 
and analyzed. Additionally, all class statistics were used simultaneously 
to classify all the selected random test fields, producing overall classi- 
fication accuracies. (See section 4.0 for examples.) 
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3.1 LANDSAT-1 Data 


The LANDSAT-1 frames over the study site were used for analysis and 
covered three distinctive phenological dates in the Southeastern United 
States. These dates covered winter (I.D. 1127-16253, November 27, 1973), 
early spring (I.D. 1217-16254, February 25, 1973), and late spring (I.D. 
1289-16254, May 8, 1973). Summer and autumn data were not available for 
analysis. These three data sets were registered and composed to form an 
11-channel (12 less 1 because of poor quality) temporal data set of a size 
roughly 500 scan lines and 600 pixels per line. The ordering and spectral 
coverages of these temporal channels are listed in table II. A black-and- 
white rendition of channel 6 (February band 2) is shown in figure 3. 

Because of the small spatial resolution of LANDSAT-1 MSS data (at 
80x80 meters per pixel after registration), the sizes of training fields 
and test fields were constrained to be no more than 5x5 pixels. Hardwood 
features usually called for narrower and smaller fields. In this study, 
the number of training fields and test fields had been chosen roughly 
proportional to their occurrence in the study site. These fields consti- 
tuted roughly 1 percent of the study site. 


3.2 Separability Measures: Divergence and PCC 


TwO mathematical measures were used to quantify the spectral separa- 
bility between timber types: (1) the divergence measure, J and (2) the 

probability of correct classification, PCC. (Marill and Green, 1963; 
Anderson, 1958; and Chang, 1971) 

The divergence measure is an approximate measure for separability, 
while the PCC measure is truly the separability measure. However, PCC is 
difficult, if not impossible to calculate except by straightforward estima 
tion via computer classification which provides classification accuracy 
measures. On the other hand, the divergence measure can be algorithmical 1 
computed based on the usual statistical assumption of normality. In this 
study, the divergence measure and the classification accuracy measure were 
jointly calculated where an increase in values of either measure was con- 
strued as an increase in separability 

The definitions of the divergence and the pairwise correct 

classification accuracy PCCA(C^,C ? ) between the two statistical classes 
and are as follows: 

J(Ci,C 2 ) - - S 2 ] [Sj - 1 - sy 1 ]] 

* 1 >[ M i • m 2] T pi' 1 * V 1 ] [ M 1 - m 2] 
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where M^, are the mean vector and covariance matrix of , and 
M 2 , Sj are those for C 2 , tr stands for the trace operation on matrices, 
and T stands for the transpose operation on matrices. 

PCCA(C 1 ,C 2 ) - ^sPCCCCj) v *sPCC(C 2 ) 

where PCC(Cj) and PCCCC 2 ) are obtained by a two-class classification 
between and C 2 , and where 

PCCfC 1 

v 1' = (number of pixels of correctly classified as from C^) 
t (number of pixels of C^) 
and PCC(C^) is defined similarly. 


When more than two classes are involved in classification, the overall 
average pairwise correct classification can be calculated as the average 
of the n(n-l)/2 pairwise measures, obtained from n(n-l)/2 possible pairs 
from the n classes (n>2). 


4.0 ANALYSIS OF DATA PROCESSING RESULTS 


All the data processing was performed at the Johnson Space Center, 
NASA, and on the Earth Resources Interactive Processing System (ERIPS) 
which is an interactive computer system developed at the center for remote 
sensing applications. This system has the capability of training field 
classification analysis as described in section 3.0. However, in this 
application on ERIPS, the transformed divergence J' was used instead of 
the divergence J defined in section 3.2, where J’=999 (l-exp(-J/16)) . 

J' and J are equivalent (Swain, 1973); any conclusion drawn from J' 
computations applies to J computations and vice versa. Without further 
complication, the following sections will abuse the notation, using J to 
denote the transformed divergence and using divergence to mean transformed 
divergence . 


4.1 Spectral Signature Plots 


Before presenting the quantitative separability results in the next 
three sections, a most effective qualitative analysis could be made by 
plotting the spectral signatures of all the forest features. A spectral 
signature plot means the graph of the statistical mean data values 
versus spectral channel. The mean values ’ere obtained from analyzing 
training fields. 
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Figure 4 is such a plot for all the 10 level V (the most detailed 
level) forest features. Four groups of features seem readily distinguish- 
able and are thus presented in the figure: (1) regeneration, nonstocked; 

(2) regeneration, seedling and sapling; (3) hardwood, immature sawtimber; 
and (4) softwood and mixed softwood -hardwood, comprising loblolly and 
shortleaf, sawtimber and poletimber, mature and immature (table I). 

In qualitative terms, then, only level II forest features can be dis- 
tinguished from one another, except for the mixed softwood-hardwood feature 
in this level. The more detailed detection levels in levels III, IV, and 
V seem too much to ask of the LANDSAT-1 MSS sensor. Additionally, temporal 
channels 6 and 11 show widest spread of data valu~s in the above mentioned 
four groups of signatures, indicating that these two channels would likely 
be the best two channels for discrimination of level II forest features 
(except for the mixed feature). Determination of the truly best channels 
involves consideration of the spread of data values about the mean value of 
all the features. Figure 4 does not indicate this kind of statistical 
variation, and only an analysis as in section 4.3 will give the most defini- 
tive answers. 

These signatures were analyzed from a total of 25 training fields of 
the 10 level V forest features. 


4.2 Pairwise Separability of Forest Features 


The 25 level V training fields selected earlier were aggregated into 
10 level V training classes, seven level IV training classes, five level III 
training classes, and four level II training classes. Each "training class" 
at any hierarchy level is taken as representative of the forest feature 
regarding spectral characteristics. Pairs of classes at all levels were 
classified one pair at one time, to obtain the pairwise correct classifica- 
tion accuracies (PCCA) . 

Figure 5 shows in bar-charts the PCCA versus feature pairs for level 
II, III, and the special level of general age classes displayed in sets of 
four bars: (1) best two channels of temporal data. (2) best two channels 

of November data, (3) best two channels of February data, and (4) best 
two channels of May data. The best channel sets were taken from results 
of section 4.3. Levels IV and V PCCA's are not presented in this paper 
but are available from Dillman and Kan (1975). PCCA between the softwood 
features of levels IV and V are generally between 50 and 60 percent, and 
those plots do not offer additional conclusions on the separability between 
those features. Plots of PCCA's for other channel set sizes are also 
available from Dillman and Kan (1975), and they show similar trends as in 
figure 5. 

Level II PCCA's are high (from 87 to 99 percent) for all feature pairs 
at the best choice of data set (i.e., seasons or combinations of seasons) 
for channel set size of 2, except that softwood can only be separated from 
mixed softwood-hardwood with less than 80 percent. The lower PCCA between 
softwood and mixed softwood-hardwood is cc.istrued to be due to the defini- 
tion of the two features (cf. section 2.2). Level III PCCA's are basically 
as high as level II PCCA's for softwood versus hardwood versus regeneration. 
Within the softwood and between softwood and mixed softwood-hardwood, PCCA's 
are as low as 61 percent. 
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The plot of figure 5(b) for the special level of general (softwood) 
age classes shows that separation between sawtimber and poletimber is poor 
(from 56 to 71 percent) where sawtimber versus regeneration or poletimber 
versus regeneration has PCCA above 95 percent. In other words, using 
LANDSAT - 1 MSS sensor, the only significantly detectable age difference is 
between an established forest versus a young (or denuded) forest, i.e., 
regeneration. 

The general trend of best season or combination of seasons is also 
discernible in figure 5; however, such conclusions are deferred to section 
4.3, where additional divergence calculations and analysis are made. 


4.3 Best Spectral Bands and Seasons 


Using the training class statistics, the best channels from the 
temporal data set and the three individual season data sets were found for 
channel set sizes of 1, 2, 3, and 4. The prioritizing of channels was 
performed by ordering the magnitude of average divergence values between 
the training classes and was done for all levels of hierarchy. Table III 
tabulates that information. 

From table III, it is apparent that the temporal data set (winter, 
early and late spring) offers the biggest overall average separability 
between all forest features at all hierarchy levels and at all channel set 
sizes (from 1 to 4 channels). The rating for individual seasons shows that 
February (phenological early spring) and May (phenological late spring) 
are better months than November (phenologically winter) for remotely sensing 
forest features in the Southeastern U.S. In particular, the red and near 
infrared channels of February and May data are good channels — channels 6 
and 10 (sometimes 11 instead of 10) of the temporal data set. 

Also, the average separability is shown to be lower at higher hierarchy 
levels, i.e., at higher details. 


4.4 Classification Accuracy Versus Channel Set Size 


Overall classification accuracies were also obtained when all classes 
were classified at the same time (in contrast with the pairwise cross-class 
classification described in section 4.2). Accuracies were obtained for all 
hierarchy levels, data sets, training fields, and a total of 24 test fields 
which were randomly distributed in the study site (Hillman and Kan, 1975). 
Only the test field overall classification accuracies are displayed in 
figure 6 versus the channel set sizes of 1, 2, 3, and 4 for the temporal 
and February data sets. Curves for May and November data sets are similar 
and lower than the February curves. 

Increase in channel set size normally improves classification accuracy 
but it was found that the two best channels performed almost as well as 
three or four best channels. Thus the temporal data set proved superior 
to the other season data sets in classification accuracy and separability. 
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Apart from the above quantitative results, the entire study site was 
also classified and displayed in figure 7 for visual comparison with the 
unclassified MSS imagery of figure 3. The classification was performed on 
the temporal data set, using the best four channels (6, 8, 10, and 11) 
with the four level II forest features. 


5.0 CONCLUSION 


This computer-aided study has investigated quantitatively the spectral 
separability of timber types in Sam Houston National Forest of Texas, which 
is a typical forest in Southeastern United States, using LANDSAT-1 multi- 
spectral scanner data. Five hierarchy levels of mapping detail plus one 
level of general (softwood) age class were studied of three sets of data at 
winter, early and late spring. Also the temporal composite of those three 
data sets was studied. Seven conclusions are summarized as follows: 

• The LANDSAT multispectral scanner sensor could be effectively used 
to separate the forest features of softwood, hardwood, and regenera- 
tion. Pairwise correct classification of training sets ranges from 
87 to 99 percent and average correct classification for test fields 
ranges from 70 to 79 percent. 

• The only significantly detectable age difference was between an 
established forest versus a young (or denuded) forest, i.e., 
regeneration. This conclusion was drawn from experience on soft- 
wood forests. 

• The red (band 2: 0.6-0. 7 micron) and one near infrared channel 

(band 3: 0.7-0. 8 micron or band 4: 0. 8-1.1 microns) of any of 

the three seasons (winter, early and late spring) would be better 
for discrimination. 

• Phenological early and late spring could be equally good seasons 
for discrimination and would be better than the winter season. 
(Summer and autumn data were not available for analysis.) 

• A temporal analysis using early and late spring LANDSAT data could 
improve classification accuracy up to 11 percent over single-season 
analysis . 

• Analysis using the two best channels would perform almost as well 
as the four best channels for s ingle - season or temporal analysis; 
hence, a two-channel analysis could be more cost-effective. 

• It would be difficult to discriminate the forest features of mixed 

so ftwood - hardwood from the softwood feature by virtue of the 
definition of the two features: (1) softwood stand contains more 

than 50 percent softwood and (2) a mixed softwood-hardwood stand 
contains between 25 and 50 percent softwood. 

It is observed that the capability of LANDSAT-1 is limited by many 
factors, including altitude, sensor design, and spatial data resolution. 

The last factor, in particular, will influence the level of mapping detail. 

A comprehensive study on forest classification and modeling has been 
reported by Kan et al. (1975). 
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TABLE II 
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Figure 1, - Location of Sam Houston National Forest (shaded area) 
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SCHEMATIC FLOW DIAGRAM OF 
TIMBER TYPE SEPARABILITY ANALYSIS 


LANDSAT-1 
MSS DATA 
(NOV, FEB, MAY) 



TEST FIELD 
LOCATION 


Figure 2 
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RAVEN DISTRICT OF SAM HOUSTON NATIONAL FOREST, 
LANDSAT - 1 IMAGE BAND 2, .6 - .7 pm, 
FEBRUARY 25, 1973 

(CONROE UNIT, STUDY SITE OUTLINED) 



Figure 3 


<_n 

ro 


NASA $ 75 10968 

LANDSAT, SPECTRAL SIGNATURE MEAN 
PLOT FOR LEVEL Y TRAINING FIELDS 



Figure 4 
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AVERAGE PAIRWISE CORRECT CLASSIFICATION 
ACCURACIES OF SAM HOUSTON NATIONAL FOREST 
TRAINING FIELDS FOR THREE LEVELS OF HIERARCHY, 
US NG BEST 2 CHANNELS FOR 4 DATA SETS: fa) LEVEL H 
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Figure 5(a) 
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AVERAGE PAIRWISE CORRECT CLASSIFICATION 
ACCURACIES OF SAM HOUSTON NATIONAL FOREST 
TRAINING FIELDS FOR THREE LEVELS OF HIERARCHY, 
USING BEST 2 CHANNELS FOR 4 DATA SETS: 

(b) LEVEL OF GENERAL AGE CLASSES 
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Figure 5{b) 
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AVERAGE PAIRWISE CORRECT CLASSIFICATION 
ACCURACIES OF SAM HOUSTON NATIONAL FOREST 

TRAINING FIELDS FOR THREE LEVELS OF HIERARCHY, 
USING BEST 2 CHANNELS FOR 4 DATA SETS: (c) LEVEL III 
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Figure 5(c) 
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AVERAGE CORRECT CLASSIFICATION ACCURACIES 
FOR SAM HOUSTON NATIONAL FOREST TEST FIELDS 
USING LANDSAT CHANNEL SETS OF SIZE 4, 3, 2, 1 
FOR SELECTED LEVELS OF HIERARCHY AND DATES 


AVERAGE 
CORRECT 
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ACCURACIES OF 
TEST FIELDS (7.) 



Figure 6 
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LEVEL H CLASSIFICATION OF CONROE STUDY SITE 
USING FOUR CHANNELS 
6, 8, 10 r 11 TEMPORAL DATA 
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Figure 7 



